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The aqueous solubility of 14 organic solutes has been calculated from their octanol/water partition coefficient 
and from their solute activity coefficient in octanol at infinite dilution. The solute aclivity coefficients were 
calculated from the Flory-Huggins and Hildebrand-Scalchard (FH-HS) equations and were found to be in good 
agreement with the activity coefficients determined from GC speeific retention volume measurements. The cal- 
culated solubilities were in good agreement with the experimental solubilities. 
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1. Introduction 

The octanol/water partition coefficient, K o/w , and the aqueous 
solubility, C 1 , of organie compounds have been used widely 
to assess and prediet the fate of toxic substances in the marine 
environment. Both of these quantities have been successfully 
correlated with environmental partitioning phenomenon such 
as bioconcentration and sorption. These two quantities are 
thermodynamically related via eq (1) [l]. 1 
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where y|£ and yjj, are the activity coefficients at infinite di- 
lution based on volume fraction for the solute in water and 
in n-octanol, respectively. The quantity, 7J£, is generally 
defined as the reciprocal of the volume fraction of the solute 

1 
in water, i.e. — and is expressed as the reciprocal of the 

<P 
solute concentration in water (C') times the molar volume of 
the solute (V), IfC'V. 

Hansch, Quinlan, and Lawrence [2] have developed a 
method for estimating /£„,„ based on the additivity principle 
of thermodynamic properties. On this basis group contribu- 
tions to log K^ or II values are defined by eq (2). 



II, = log K^Jx) - \ogKjJJi) 



(2) 



where K o/w (x) is the partition coefficient of a derivative and 
Ko/JJi) is that of the parent molecule. Despite evidenee that 
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II values are often not additive, usually for steric reasons but 
sometimes through electronic and hydrogen bonding effects, 
there are many series of compounds in which II appears to 
be invariant. 

Equations (1) and (2) suggest that the aqueous solubility, 
C', may be calculated from the molecular structure of the 
solute provided y* can be estimated with some degree of 
accuracy. There are methods [3] available for calculating 
yjj, from the solute physical properties. Because of the sim- 
plicity of the calculations, the Flory-Huggins and Hildebrand- 
Scatchard (FH-HS) method was chosen to evaluate yjj,. 

According to solution theories |4], the solute activity coef- 
ficient (y x ) based on mole fraction consists of an athermal 
term and a thermal term. At infinite dilution, in terms of the 
Flory-Huggins equation for the athermal term and the Hil- 
debrand-Scatchatd equation for the thermal term, this may 
be expressed [5] as 



Iny™ = In — h 
r 
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where r is the ratio of the molar volume of the solveut V" to 

that of the solute V, I — 1 , and the 5 and 8' are, respectively, 

the solute and solvent solubility parameters. These are de- 
fined [6] as 



S = 
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where &E V and AH V are, respectively, the molar energy and 
enthalpy of vaporization. At infinite dilution the volume frac- 
tion activity coefficient (y^) may be expressed in terms of 
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mole fraction activity coefficient (-y*) by the following equa- 
tion [ 5, 7]. 



■y* = -v, 



(5) 



Now combining eq (3) and (5), the solute activity coefficient 
"V J may be expressed as 



^=(1-; + ^(S'-S) 2 



(6) 



In this paper we are concerned with estimating the solute 
activity coefficients in water, "yjjj, using eq 1, from the esti- 
mated values of -yjj> an d me experimental values of K o/w . 
Activity coefficients, "y£, in octanol for a number of solutes 
have also been calculated from the Flory-Huggins and Hil- 
debrand-Scatchard equations and compared with the exper- 
imental values obtained by the gas chromatographic (GC) 
method. 

2. Experimental Procedure 

The stationary phase, n-octanol, used in this study was 
obtained from the Aldrich Chemical Company. The purity 
was checked by GC and was found to be 99 + mole %. The 
stationary phase was coated onto the support material, Chro- 
mosorb W-HP, 100/120 mesh, and the weight percent of 
coating was determined by an ashing method [8]. A stainless 
steel column (1/8 in OD) was then packed with a known 
amount of the coated support and was connected to a Hewlett- 
Packard 5830A 2 gas chromatograph equipped with a flame 
ionization detector (FID) for the determination of the solute 
retention time. The temperature of the column was controlled 
by circulating water through a copper tubing jacket around 
the column and by a Haake Model FK temperature regulator 
which regulated the water temperature to 25.0 ± 0.05 °C. 

The column inlet pressure, measured with a precalibrated 
pressure gauge (range 0—15 PSI), was kept constant during 
a run by regulating with a precision valve. In order to measure 
the carrier gas flow rate the column was disconnected from 
the FID, just before and after the experiment, and a soap 
bubble flowmeter was connected to the column outlet. The 
carrier gas was presaturated with n-octanol in order to reduce 
bleeding of the stationary phase. 

Since the retention times decreased as the experiment pro- 
gressed due to column loss of the rt-octanol, toluene was used 
as a reference retention time standard and injected with each 
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specify the experimental procedure, [u no case docs sued identification imply rec- 
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solute. The measured retention times were then corrected for 
bleeding using the reference solute retention time. 

The solute specific retention volumes (K°) were calculated 
by using the following equation [9]. 



VI = 



t R -F P a - P„ 273.2 3 (P/PJ 2 - 1 
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where t R is the retention time (min), F is the carrier gas flow 
rate (ml/min), P t and P a are the column inlet and atmospheric 
pressure (in Torr), P w is the saturated water vapor pressure 
(in Torr) at room temperature T ( K), and W is the weight (in 
gram) of stationary phase packed in the column. The reported 
specific retention volumes listed in column 1 of table 1 are 
an average of 6 measurements. 

3. Results and discussion 

The solute activity coefficients at infinite dilution (7^) were 
determined from the following equation [8]. 



In-y* = In 
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where M' is the molecular weight of stationary phase, n- 
octanol, (g/mole), P° and B are, respectively, the vapor pres- 
sure (in Torr) and the gas phase second virial coefficient (ml/ 
mole) of the pure solute at temperature T^K), and R is the 
gas constant. The second virial coefficients were calculated 
from the McGlashan and Potter equation [10] using the con- 
stants tabulated by Dreisbach [11, 12, 13] and Timmermans 
[14]. These sources were also used for Antoine constants to 
calculate values of P° and for constants for the law of recti- 
linear diameters to compute values of molar volumes. The 
solute mole fraction and volume fraction activity coefficients 
listed in table 1 were calculated using equations 8 and 5. It 
is estimated that the probable uncertainty in the activity coef- 
ficients listed in table 1 is about 2 %. 

The solute solubility parameters, 8, listed in table 1 were 
calculated from their molar enthalpies of vaporization listed 
in the Dreisbach compilations [11, 12, 13]; whereas the 
solvent solubility parameter, S' was treated as an adjustable 
parameter. In order to fit the experimental data of the systems 
studied, two different solubility parameters (8') were used, 
one for aromatic solutes (8' = 6.71) and another for aliphatic 
solutes (8' = 9.57). The solute activity coefficients calcu- 
lated using equation 6 are listed in column 5 of table 1. The 
agreement between the two sets of data (calculated and ex- 
perimental activity coefficients) is in general good. The av- 
erage deviation is 8.7%; however, in a few systems it is as 
high as 18%. These results suggest that the activity coefficient 
o( a system may be calculated from the knowledge of their 
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Table 1. Specific Retention Volumes (P'°) and Solute Activity Coefficients in n-Oclanol at 25.0 °C. 





^ 




tS 


B- 


T% 


TZ(e* v l) - -yKcalc) ^ ioq 


Solute 


(ml/g) 


% 


(expt) 


(cal/cm 3 )" 2 


(calc) 


7i(expl) 


re-Pentane 


100.6 


2.62 


3.59 


7.17 


4.04 


-12:5 


K-Hexane 


312.5 


2.81 


3.39 


7.33 


3.61 


-6.5 


n-Heptane 


943.0 


3.06 


3.29 


7.47 


3.21 


2.4 


re-Octane 


2793. 


3.36 


3.25 


7.55 


2.99 


8.0 


1-Hexene 


280.5 


2.54 


3.19 


7.38 


3.40 


6.6 


1-Heptene 


— 


— 


— 


7.53 


3.01 


— 


1-Octene 


2518. 


3.00 


3.01 


7.61 


2.79 


7.3 


1-Nonene 


7589. 


3.24 


2.95 


7.73 


2.45 


16.9 


benzene 


669.9 


2.07 


3.67 


9.20 


3.94 


-7.4 


Toluene 


2120. 


2.18 


3.23 


8.91 


3.32 


-2.8 


Ethylbenzene 


5542. 


2.49 


3.20 


8.79 


3.07 


4.1 


o-Xylene 


8317. 


2.38 


3.12 


8.99 


3.67 


-17.6 


m-Xylene 


6763. 


2.34 


3.00 


8.82 


3.15 


-5.0 


"-Propylbenzeiie 


13854. 


2.76 


3.12 


8.64 


2.71 


13.1 


Isopropyl- 














benzene 


10623. 


2.66 


3.00 


8.61 


2.64 


12.0 



1 calorie = 4. 184 Joules 



physical properties using the Flory-Huggins and the Hilde- 
brand-Scatchard equations. 

The quantity "yjj, in eq 1 is defined as the activity coefficient 
of a solute measured in octanol saturated with water. How- 
ever, Purnell [15] has shown that 7JJ, for a solute in a binary 

mixture may be expressed as 



Table 2. Activity Coefficients and Octanol/Water Partition Coefficients at 
25.0 X. 









+ 






(9) 



where 7^ and 7^' are the solute activity coefficients at 
infinite dilution in pure octanol and in pure water, respec- 
tively, and <t and <$> w are their volume fractions. It has beeu 
reported [16] that at equilibrium the volume fraction of water 
in octanol is 0.0414. The value of 7$ for the solutes used in 
this study vary between 380 and 63,000 while those of 7^ 
vary between 3.00 and 3.60. Thus the errors involved in 
assuming 7^ = 7^ and 7^ = 7JJ are minimal. 

The solute activity coefficients in water, 7^, listed in col- 
umn 3 of table 2, were calculated using the experimental log 
K^ and calculated log 7$ from the Flory-Huggins and Hil- 
debrand-Scatchard equation. Considering that the average 
standard deviation in the experimental log K a/w is 0.04 and 
log 7^ is 0.03, the agreement between the experimental log 
7^ (listed in the last column) and the calculated log 7^ is very 
encouraging. 

In summary this study suggests that the solute activity 
coefficient 7J, may be calculated accurately using the Flory- 
Huggins and Hildebrand-Seatchard equations, and the sol- 
ubility of organic eompounds in water may be predicted suc- 
cessfully from the knowledge of their octanol/water partition 





'og-y* 


log JW17) 


log 11 


Experimental 


Solute 


calculated 


Experimental 


calculated 


(17) 


n-Pentane 


0.61 


3.62 


4.23 


4.19 


n-Hexane 


0.56 


4.11 


4.67 


4.73 


n-Heptane 


0.51 


4.66 


5.17 


5.27 


n-Octane 


0.48 


5.18 


5.66 


5.80 


1-Htxene 


0.53 


3.39 


3.90 


3.99 


1-Heplene 


0.48 


3.99 


4.47 


4.58 


1-Oetene 


0.45 


4.57 


5.02 


5.24 


1-Nonene 


0.39 


5.15 


5.54 


5.82 


Benzene 


0.60 


2. 13 m 


2.73 


2.59 12 ' 


Toluene 


0.52 


2.65 


3.17 


3.17 


Ethyl- 










benzene 


0.49 


3.13 


3.62 


3.69 


o-Xylene 


0.56 


3.13 


3.69 


3.60 


m-Xylene 


0.50 


3.20 


3.70 


3.73 


re-Propyl- 










benzene 


0.43 


3.69 


4.12 


4.22 


Isopropyl- 










benzene 


0.42 


3.90 


4.32 


4.09 



coefficients and the solute activity coefficient calculated from 
the FH-HS equation. 
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